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ABSTRACT 

Remotely sensed ocean color products require highly accurate top-of-atmosphere (TOA) radiances, on the order 
of 0.5% or better. Due to incidents both prelaunch and on-orbit, meeting this requirement has been a consistent 
problem for the MODIS instrument on the Terra satellite, especially in the later part of the mission. The NASA 
Ocean Biology Processing Group (OBPG) has developed an approach to correct the TOA radiances of MODIS 
Terra using spatially and temporally averaged ocean color products from other ocean color sensors (such as the 
SeaWiFS instrument on Orbview-2 or the MODIS instrument on the Aqua satellite). The latest results suggest 
that for MODIS Terra, both linear polarization parameters of the Mueller matrix are temporally evolving. A 
change to the functional form of the scan angle dependence improved the quality of the derived coefficients. 
Additionally, this paper demonstrates that simultaneously retrieving polarization and gain parameters improves 
the gain retrieval (versus retrieving the gain parameter only). 
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1. INTRODUCTION 

Two copies of the Moderate- Resolution Imaging Spectroradiometer (MODIS) 1 are currently in operation, provid- 
ing top-of-atmosphere (TOA) radiances from 412nm to 14200nm. The first unit was launched in December 1999 
on NASA’s Earth Observing System (EOS) Terra satellite, the second on the Aqua satellite in May 2002. MODIS 
has 36 spectral bands on four different focal planes. The Ocean Biology Processing Group (OBPG) at NASA 
uses both instruments to produce ocean color data products 2 and provides the data (approximately 23 Gigabyte 
per day for each instrument) to users around the world from the following website: http:oceancolor.gsfc.nasa.gov. 

Each MODIS ocean color band has 10 independent detectors, and each of these detectors needs to be calibrated 
and characterized separately. Both MODIS instruments are calibrated using on-board calibrators 3 and lunar 
irradiances. 4 For MODIS Terra, these calibrations have not been sufficient to produce high quality ocean color 
products. 5 The OBPG has developed a cross calibration method to improve the MODIS Terra characterization, 
using water-leaving radiances from another sensor as truth fields. 6, 7 This paper presents an update of the derived 
corrections relative to the MODIS Aqua instrument. Note that for the operational data processing, MODIS Terra 
ocean color products use SeaWiFS 8 for the early part of the mission, and MODIS Aqua for the later part of the 
mission (SeaWiFS data is only available until 2010). The focus of this paper is the presentation of the temporal 
evolution of the MODIS Terra polarization sensitivity, which did not change significantly before the availability 
of MODIS Aqua data, so it was decided to use only MODIS Aqua data as truth data for this paper. 

2. STANDARD MODIS CALIBRATION AND CHARACTERIZATION METHODS 

The standard MODIS calibration equation 9 uses the calibration coefficient mi to describe changes in the radio- 
metric sensitivity of the instrument. The mi’s are used to derive earth scene reflectance pev from the measured 
counts of the earth scene ( driEV ) with: 

Pev ■ cos(9 E v) = mi ■ dn* EV ■ d 2 Earth _ Sun (1) 
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where pev is observed earth-viewing reflectance at a MODIS spectral band, dn* EV is the temperature corrected 
measured counts, and dEarth-Sun is the distance between earth and sun, and 6ev is the solar zenith angle. The 
scan angle dependence is corrected by dividing the measured counts by the RVS (response- versus-scan). 

The calibration source for determining mi on-orbit is the solar diffuser. The solar diffuser is viewed at 
an angle of incidence on the scan mirror of 50.3°. The calibration source for determining the RVS are lunar 
measurements through the space view port. For some bands only, measurements of selected earth view targets 
(several desert areas) are used for gain trending, see Sun et al . 10 

3. VICARIOUS CALIBRATION METHOD 

The standard polarization correction equation for an uncalibrated instrument 11 is 

L m = MnL t + Mi 2 (QtCOs2a + U t sin2a) + Mi 3 (—Q t sin2a + U t cos2a ) + M u V t (2) 

where (Lt, Qt , Ut, V t ) is the Stokes vector at the TOA, L m is the measured radiance, Mn to M 14 are the elements 
of the first line of the Mueller Matrix, and a is a rotation angle to adjust for different reference frames of the 
modeled atmospheric polarization and the instrument. Since Vt, the circular polarization component, is very close 
to zero at the TOA, M 13 is irrelevant. The polarization sensitivity parameters M 12 and M 13 were determined 
prelaunch at scan angles from —45° to +45°, for each band, mirror side, and detector . 12 The small variations of 
these parameters with detector were considered suspect and not applied in the ocean color processing . 13 Note 
that Mn = 1-0 means that eq. 1 (including the RVS corrections) was adequate, and M 32 = M 13 = 0 means that 
there is no polarization sensitivity. 

The cross calibration method has been described by Kwiatkowska et al ., 6 so here we provide only a brief 
summary. The correction parameters Mu, M 12 , and M 13 were calculated relative to the calibration (version 
6.1.17.17„OC) provided by the MODIS Characterization Support Team (MCST). For a given day, the level 3 
water-leaving radiances from MODIS Aqua are used to predict the top-of-atmosphere (TOA) radiances as seen 
by MODIS on that day, using the vicarious calibration technique described in Franz et al., 2007, 14 and the 
atmospheric correction model consistent with the R2010.0 reprocessing. Thus, all relevant components of the 
Stokes vector ( L t , Q t , Ut ) can be modeled. This allows the estimation of the radiometric calibration parameters 
Mn as well as the polarization sensitivity parameters M 12 and M 13 by minimizing differences between the 
left-hand and right-hand side of eq .2 over many observations. 

The parameters Mu, Mi 2 , and M 13 are derived simultaneously by constructing a matrix of thousands of 
instances of eq. 2, each line corresponding to a matchup of one MODIS Terra measurement ( L m ) with one 
predicted TOA Stokes vector. 

In the correction process, Mu is applied first, i.e. Mi 2 and M 13 need to be applied relative to a calibrated 
instrument. The polarization correction equation for a calibrated instrument 11 is thus 

L m = L t + m\ 2 {QtCos2a + Utsin2a) + m\ 3 (—Qtsin2a + UtCos2a) + muVt (3) 


with 


mi 2 

mi 3 


Mi 2 

Wi 

M13 

Mn 


( 4 ) 

( 5 ) 


The modeled TOA radiances are compared to the radiances measured by MODIS for every scan angle, mirror 
side, and detector. This means that the instrument characterization parameters (Mn, mi 2 , m. 13 ) can be derived 
as a function of scan angle, mirror side, and detector. To reduce noise, the scan angle dependence is modeled by 
a fourth order polynomial for Mn, and as a quadratic function for mi 2 and mi 3 . 

The method assumes that the near-infrared (NIR) bands do not need a correction. Fig. 1 shows that the 
globally averaged aerosol optical thickness (AOT) ocean color products have very similar temporal trends for 
MODIS Terra and MODIS Aqua, which suggests that the MCST calibration of the MODIS NIR bands indeed 
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Figure 1. Globally averaged MODIS AOT (left) and angstrom (right), for MODIS Aqua (solid line) and MODIS Terra 
(dashed line). The methodology to create these plots is explained in Franz et al. 5 


does not need a correction. The Angstrom trends for Terra and Aqua MODIS are generally in good agreement, 
but with some bias likely associated with the vicarious calibration. The fact that both the AOT and Angstrom 
trends remain in good agreement after 2005, when the estimated polarization changes in MODIS Terra become 
significant, suggests that the polarization change detected in the visible bands is not an artifact of the NIR 
calibration. 

The retrieval of the instrument characterization parameters is repeated for one day in every month of the 
mission. This results in a time series for the instrument characterization parameters. The parameters are 
retrieved as a function of scan angle for each day. Examples for Mu, TO12 and TO13 of band 8 (412nm) are shown 
below in Figs. 2, 4 and 5, respectively, to give an estimate of the noise in the measurements. 

In Fig. 3, Mu were retrieved while keeping mi2 and m.13 fixed at their prelaunch values. It can be seen 
that at the high scan angles, not including the polarization parameters in the retrieval more than doubles the 
amplitude of the seasonal oscillation for the Mu- Furthermore, the long term trends have changed significantly 
(e.g. for frame 1250 at the end of mission, the Mu in Fig. 2 are at about 1.1, versus only about 1.05 for Fig. 3. 
A difference of 5% in the gain calibration is highly significant, not only for ocean color applications. 

The MOST approach does not take into account a changing polarization sensitivity. Neither do the results 
shown in Fig. 3. As expected, the Mu in Fig. 3 are closer to 1.0 than the Mu in Fig. 2. The large deviations of 
the Mu from 1.0 are therefore likely a result of the fact that the MOST calibration approach does not take into 
account polarization sensitivity, and not an indication that there is systematic problem with their approach for 
unpolarized radiances. 

A. Lyaspustin (private communication) reported a potential overcorrection at the beginning of the scan when 
applying the operational OBPG correction coefficients to MODIS Terra land products (see also Lyapustin 15 ). 
We discovered that the reason was that the scan angle dependence for the 77112 coefficients was modelled as linear, 
but a quadratic function should have been used. The large increase of polarization sensitivity at the end of scan 
produced a large slope of the linear function, which lead to mi2 values at the beginning of the scan that were 
too low, see Fig. 2 in Meister et al. 7 (the mi2 decrease to about -0.1 after 2007 for band 8 for the beginning of 
the scan). The data presented here use a quadratic function to fit the scan angle dependence of mi2 and TO13. 
As can be seen in Fig. 4, the large increase in mi 2 at the end of scan (frame 1250) no longer produces a decrease 
of the TO12 at the beginning of scan (frame 22). In effect, there is now a small increase of TO12 for frame 22 after 
2007. 
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Figure 2. Mn for mirror sides 1 (solid line) and 2 (dashed line), detector 1-10 (black to red) of MODIS Terra band 8 
(412nm)as a function of time at four different scan angles (see plot title). All parameteres (Mn,mi2 and 77113) were 
retrieved simultaneously. Note the different ordinate (y-axis) for each plot. 
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Figure 3. Mn for mirror sides 1 (solid line) and 2 (dashed line), detectors 1-10 (black to red) of MODIS Terra as a 
function of time at four different scan angles (see plot title). Only Mn was retrieved, m \2 and TO13 were fixed at their 
prelaunch values to demonstrate that retrieving both gain and polarization paramters simultaneously improves the gain 
retrievals. 
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Figure 4. 11112 for mirror sides 1 (solid line) and 2 (dashed line), detector 1-10 (black to red) of MODIS Terra band 8 
(412nm) as a function of time at four different scan angles (see plot title). All parameteres (Mn,mi 2 and m\z) were 
retrieved simultaneously. Note the different ordinate (y-axis) for each plot. 
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Figure 5. mi 3 for mirror sides 1 (solid line) and 2 (dashed line), detector 1-10 (black to red) of MODIS Terra band 8 
(412nm) as a function of time at four different scan angles (see plot title). All parameteres (Mn,mi 2 and 771 , 13 ) were 
retrieved simultaneously. Note the different ordinate (y-axis) for each plot. 


4. RESULTS AFTER TEMPORAL SMOOTHING 


It can be seen in Figs. 2, 4 and 5 that there is significant noise in the retrievals. To reduce noise, a running 1 year 
median is calculated for 15 different scan angles (for equally spaced scan angle intervals), and the polynomial as 
a function of scan angle is recalculated for the median values for each date. To further reduce noise, a temporal 
smoothing filter is applied (± 2 months smoothing). This method was consistenly applied to all bands and all 
three parameters (Mu, m \2 and mis). Previously, (see e.g. Meister et al. 7 ) the m\s were kept constant at their 
prelaunch values. The additional amount of data available since the latest analysis gave us the confidence to 
trend the TO13 as well. 

An overview of the results for mirror side 1, detector 5, is given in Figs. 6 to 8. The shortest wavelengths 
show the largest values, and the end of scan varies much more than the beginning of the scan for most bands. 
These effects are expected, because the shortest wavelengths have always been more challenging to calibrate for 
the MODIS instruments, and the end of scan has the largest angle of incidence on the scan mirror, which favors 
a stronger polarization sensitivity. For band 10 (488nm), the largest correction for Mu is actually in the center 
of the scan. Note that the results for 748nm are just a consistency check: the NIR bands are assumed to be 
perfectly calibrated by our approach, i.e. the Mu should be close to 1.0 and the mi2 and 77113 should be close 
to their prelaunch values; for band 16 (869nm, not shown here) the results are even better than for band 15 
(748nm). 

For almost all bands, the magnitude of the parameters increases with time. This is not surprising, the more 
the instruments (both MODIS Aqua and Terra) age, the larger correction coefficients we expect. 

The results for all detectors and both mirror sides at the end of scan are shown in Figs. 9 to 11. Fig. 9 shows 
the result for Mu . It can be seen that especially for the shorter wavelengths, there is a clear mirror side difference. 
For the longer wavelengths (555nm to 678nm), the mirror sides agree quite well, but one detector (detector 1, 
black line) is an outlier. In general, shorter wavelengths show larger values (up to almost 1.2 for 412nm and 
443nm). The red MODIS ocean bands (667nm and 678nm) show very good agreement between MODIS Aqua 
and Terra (within ±1%), but the red MODIS Terra land band (645nm) shows an increase of Mu of about 5% 
(there is a similar increase at 859nm). The relatively large difference between MODIS Terra and MODIS Aqua 
at 645nm (and 859nm) is rather surprising, it does not appear in independent trend analysis by MOST of those 
bands at high radiances, but was confirmed by MOST for ocean radiances (private communication, Amit Angal, 
MOST). This could mean that the linearity of those bands is changing with time. 

Fig. 10 shows the result for 77112. It can be seen that the shortest wavelength bands had a very large increase 
in polarization sensitivity over time after 2005. Mirror side 2 starts to increase earlier than mirror side 1. At 
412nm, 77112 of mirror side 2 reaches a maximum of about 0.4 and then decreases. No other band shows a similar 
behavior. The 645nm band has almost no change in polarization sensitivity over time, all other bands change. 

Fig. 11 shows that not only band 8 (412nm) shows a decrease of 77113, but also some longer wavelengths, 
especially bands 13 and 14 (667nm and 678nm). The amount of the decrease is up to 0.02. 

Note that some of the plots in Figs. 9 to 11 show short term variations that may or may not be valid changes 
of the MODIS Terra characteristics. For a reprocessing of the MODIS Terra data, the OBPG may decide to 
apply further temporal smoothing before applying corrections to the MODIS Terra data. 

5. CONCLUSIONS 

Recently, there has been increased interest from other disciplines (land, atmosphere) in the MODIS Terra po- 
larization correction, and several studies( 15 ^ 17 ) have reported increased product quality using the corrections 
derived by the OBPG. This paper provides a reanalysis of the correction factors for MODIS Terra and extends 
the time period covered through 2013. Note that the OBPG expects to reprocess the truth field used in this 
analysis (the MODIS Aqua data record) in the near future. Therefore, the coefficients presented here are only 
preliminary, but they do provide insight into the expected changes. 

Extending the time period showed that the polarization correction parameter 777.13 is not constant, as was 
assumed previously. 6, ' The application of this change has not yet been tested in the correction of ocean color 
products, but it will be in the near future. 
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Figure 6. Mu for mirror side 1, detector 5 of MODIS Terra as a function of scan angle. Colors indicate time (black=2002, 
red=2013). Note the variation of the ordinate axis among the plots. 
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Figure 7. mi 2 for mirror side 1, detector 5 of MODIS Terra as a function of scan angle. Colors indicate time (black=2002, 
red=2013). Note the variation of the ordinate axis among the plots. 
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Figure 8. mi 3 for mirror side 1, detector 5 of MODIS Terra as a function of scan angle. Colors indicate time (black=2002, 
red=2013). 
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Figure 9. Mu for mirror sides 1 (solid line) and 2 (dashed line), detectors 1-10 (black to red) of MODIS Terra as a function 
of time for all wavelengths at the end of scan scan (frame 1250). 
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Figure 10. 1 U 12 for mirror sides 1 (solid line) and 2 (dashed line), detectors 1-10 (black to red) of MODIS Terra as a 
function of time for all wavelengths at the end of scan scan (frame 1250). 
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Figure 11. mi 3 for mirror sides 1 (solid line) and 2 (dashed line), detectors 1-10 (black to red) of MODIS Terra as a 
function of time for all wavelengths at the end of scan scan (frame 1250). 


The functional form of the scan angle dependence for the polarization correction parameters m \2 and mi 3 was 
changed from linear 6, 7 to quadratic. This change eliminated an out-of-family behavior of the m \2 coefficients 
of band 8 (412nm) at the beginning of the scan in the later half of the Terra mission and should reduce a 
potential overcorrection reported in the application of the polarization correction coefficients to MODIS Terra 
land products. 15 

For most bands, the polarization sensitivity increases with time. A notable exception is band 8, mirror side 
2, where the polarization correction parameter m \2 increased up to 2009 and then fell in the following years to 
less than half its maximum value. This rather surprising temporal evolution corresponds to an unusual temporal 
trend of the response versus scan for band 8 mirror side 2 (not shown here). 

A comparison of the cross calibration approach retrieving both gain corrections and polarization correction 
parameters versus retrieving the gain corrections only showed that not including the polarization correction 
parameters resulted in a large increase of the seasonal oscillation in the retrieved gain corrections for band 8 
and in significant differences in the long term trends as well. This underscores the importance of deriving and 
applying both gain corrections and polarization correction parameters simultaneously. 
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